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INTRODUCTION 


The  elastic-plastic  problems  of  thick-walled  tubes  subjected  to  mechani¬ 
cal  loadings  have  been  solved  by  many  investigators  based  on  different  theo¬ 
ries  or  methods. Very  little  work  has  been  done  on  the  elastoplastlc  solu¬ 
tions  for  thick-walled  tubes  subjected  to  thermal  loadings . Recently,  a 
new  finite-difference  approach  was  developed  for  solving  the  plane-strain 
problems  of  elastic-plastic  thick-walled  tubes  subjected  to  mechanical  load¬ 
ings  3  or  thermal  loadings.?  The  approach  has  been  extended  to  solve  the  gen¬ 
eralised  plane-strain  problem  subjected  to  mechanical  loadings. ^ 

In  the  present  report,  the  generalized  plane-strain  problems  of  elastic- 
plastic  thick-walled  tubes  subjected  to  mechanical  as  well  as  thermal  loadings 
are  considered.  The  formulation  includes  internal  pressure,  external  pres¬ 
sure,  axial  force,  steady  or  transient  thermal  loadings.  The  numerical  result 
is  reported  for  a  closed-end  tube  subjected  to  internal  pressure  and  to  tem¬ 
perature  gradient.  The  formulation  is  based  on  the  incremental  finite- 
difference  method  using  von  Mlses’  criterion,  the  Prandtl-Reuss  flow  theory 
and  the  isotropic  hardening  rule.  All  incremental  quantities  are  determined 
in  the  program  and  no  iteration  is  needed.  In  order  to  improve  the  efficiency 
of  the  program,  a  scaled  loading  approach  has  been  implemented.  The  numerical 
results  have  been  compared  with  those  of  Bland^  and  additional  results  due  to 
large  temperature  gradient  are  reported. 


References  are  listed  at  the  end  of  this  report. 
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BASIC  EQUATIONS 


Assuming  small  strain  and  no  body  forces  in  the  axisymmetric  state  of 
generalized  plane  strain,  the  radial  and  tangential  stresses,  ar  and  09,  must 
satisfy  the  equilibrium  equation, 

r(3or/3r)  -  09  -  or  ;  (1) 

and  the  corresponding  strains,  £r  and  £9,  are  given  in  terms  of  the  radial 
displacement,  u,  by 

er  -  3u/3r  ,  £9  -  u/r  .  (2) 

It  follows  that  the  strains  must  satisfy  the  equation  of  compatibility 

r(3£9/3r)  -  er  -  eg  ,  (3) 

Whereas  the  differential  equations  (1),  (2),  and  (3)  hold  throughput  the  tube 
regardless  of  the  material  properties,  the  constitution  equations  assume 
various  forms  according  to  the  adopted  form  of  yield  function,  hardening  rule, 
total  or  Incremental  theory  of  plasticity.  In  the  present  report,  the 
material  is  assumed  to  be  elastic-plastic,  obeying  the  von  Mlses*  yield 
criterion,  the  Prandtl-Reuss  flow  theory  and  the  isotropic  hardening  law.  The 
complete  stress-strain  relations  are: 

Ae*'  -  A0i'/2G+  ( 3/2) 0^' Ao/(oH' )  (4) 

Ao  >  0  for  i  ■  r,0,z 

Ae*  -  E-l(l-2v)Aonl  +  aAT  (5) 

where  E,  v,  a  are  Young’s  modulus,  Poisson's  ratio,  coefficient  of  thermal 
expansion,  respectively,  AT  is  the  temperature  Increment, 


2 


V 


v 


2G  -  E/(l+v)  , 

em  “  (er+ee+ez)/3  ,  e*'  -  c*  -  %  , 

°m  “  (or+og+oz)/3  ,  oi'  -  Oi  -  <%  , 

o  -  (l//2)l(or-ae)2  +  (ae-az)2  +  (o^Of)2]1/2  >  o0  ,  (6) 

and  0O  is  the  yield  stress  In  simple  tension  or  compression.  For  a  strain- 
hardening  material,  H'  is  the  slope  of  the  effective  stress/plastic  strain 
curve 

o-HC/deP)  .  (7) 

For  an  ldeally-plastic  material  (H*  -  0),  the  quantity  (3/2)do(dH' )  is 
replaced  by  dX,  a  positive  factor  of  proportionality.  When  o  <  o0  or  do  <  0, 
the  state  of  stress  is  elastic  and  the  second  term  in  equation  (4)  disappears. 

Following  Yamada  et  al,®  equations  (4)  and  (5)  can  be  rewritten  in  an 

incremental  form 

Aoi  -  dijAej  -  EoAT/(1-2v)  for  i,J  •  r,9,x 

and 

d1:|/2G  -  v/(l-2v)  +  -  ot'oj'/S  ,  (8) 

where 

2  1 

S  -  -  (1  +  -  H'/G)c2  ,  H'/E  -  u>/(l-u>)  ,  (9) 

3  3 

uE  is  the  slope  of  the  effective  stress-strain  curve,  and  6jj  is  the  Kronecker 
delta. 

Consider  an  open-end  or  closed-end  thick-walled  tube  of  inner  radius  a 
and  external  radius  b.  The  tube  is  subjected  to  inner  pressure  p,  external 
pressure  q,  and  force  f,  inner  temperature  Ta  and  external  temperature  Tb» 

The  boundary  conditions  for  the  problem  are 
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or(a,t)  -  -p  ,  ar(b,t)  ■  -q  , 
b 

2tt J  razdr  -  pita2  p  +  f  (10) 

a 

where  p  is  0  for  open-end  tubes  and  1  for  closed-end  tubes. 

When  the  temperature  T  is  not  varying  with  respect  to  time,  the  steady 
state  distribution  is  given  by 

T  -  Ta  +  N  log  (r/a) 

and 

N  “  (Tb-Ta)/log  (b/a)  (11) 

The  stress  solution  in  the  elastic  range  is  well-known.  The  quantities  p*, 
q*,  £*,  Ta*  or  Tb*  required  to  cause  initial  yielding  can  be  determined  by 
using  the  von  Hises'  yield  criterion. 

FINITE-DIFFERENCE  FORMULATION 

For  loading  beyond  the  elastic  limit,  an  incremental  approach  of  the 
finite-difference  formulation  is  used.  The  cross  section  of  the  tube  is 
divided  into  n  rings  with  ri-a,r2 , . . . .r^-p ,  ....rn+j-b,  where  p  is  the  radius 
of  the  elastic-plastic  Interface.  At  the  beginning  of  each  Increment  of 
loading,  the  distribution  of  temperature,  displacements,  strains,  and  stresses 
is  assumed  to  be  known  and  we  want  to  determine  Au,  Aer,  Aeq,  Aez,  Aar,  Aoq, 
Aaz  at  all  grid  points  for  the  applied  incremental  loading,  Ap,  Aq,  Af,  AT^  (i 
-  1  to  n-H).  Since  the  incremental  stresses  are  related  to  the  incremental 
strains  by  the  incremental  form  (Eq.  (8))  and  Au  •  rAeg,  there  exists  only 
three  unknowns  at  each  station  that  have  to  be  determined  for  each  Increment 
of  loading.  Accounting  for  the  fact  that  the  axial  strain  ex  is  Independent 
of  r,  the  unknown  variables  in  the  present  formulation  are  (Aeo)^,  (Aer)i,  for 
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1  -  1 ,2, . .  .n,irt-l ,  and  Aez. 

The  equation  of  equilibrium  (1)  and  the  equation  of  compatibility  (3)  are 
valid  for  both  the  elastic  and  the  plastic  regions  of  a  thick-walled  tube. 

The  finite-difference  forms  of  these  two  equations  at  i  -  l,...,n  are  given  by 
(ri+i-2ri)(Acr)i  -  (ri+1-ri)( Aoe)*  +  ri(Aar)*+1 

-  ( r i+i_r i) ( o 0-or) i  -  ri[(or)i+i  -  (or)il  (12) 

for  the  equation  of  equilibrium ,  and 

(ri+i-2ri)(Aee)i  -  (ri+i-riXACr)*  +  ri(Aee)1+1 

■  (ri+i-riXer-ee)i  “  fiKe0)i+l  "  C eo>il  (13) 

for  the  equation  of  compatibility.  With  the  aid  of  the  incremental  stress- 
strain  relations  (Eq.  (8)),  equation  (12)  can  be  rewritten  as 
Kri+l“2ri)(d12H  +  ("*i+l+ri)(d22)iKAee)i 
+  [(ri+i-2ri)(du)i  +  (-r1+1+ri)(d2i)i](Aer)1 
+  ri(di2)i+l(Ae6)i+l  +  ri(dll>i+l(Aer>i+l 
+  Kri+l"2riXd13H  +  (”ti+i+ri)(d23)  i  +  *i(d13)i+ll Aez 

-  (r1+1-ri)(O0-or)i  -  ri[(or)i+1  -  <or)i] 

+  rA  Ea(l-2v)-l(ATi+l-ATi)  .  (14) 

The  boundary  conditions  for  the  problem  are 

Acr(a,t)  -  -Ap  ,  Aor(b, t)  -  -Aq  , 
n 

n  l  tri(Aoz)i.  +  ri+i(Aoz)i+1](ri+i-ri)  -  pira2Ap  +  Af  ,  (15) 

i-1 

where  p  is  0  for  open-end  tubes  and  1  for  closed-end  tubes.  Using  the 
incremental  relations  (Eq.  (8)),  we  rewrite  equation  (15)  as 
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(di2>iUee)i  +  (dn)i(Aer)i  +  (di3)iA£z  -  -Ap  +  Eo(l-2v)_1  ATj  (16) 
(d12>n+l(Ae0>n+l  +  (dll)n+l(Aer)n+l  +  (d13WlAez  “  “M  +  Ect(l-v)-1ATn+1 

(17) 


and 


n 

I  (ri+l-riHr1[(d23)i(Ae0)1  +  (d13)1(  Ae,.)!]  +  r1+1[  <  d£3>  i-»-i  ( Aee)1+1 
i-1 

n 

+  (d13)i+l(Aer)i+i] }  +  l  (ri+1-r1)[r1(d33)1  +  ri+1(d33)1+1] Aez 

i-1 

n 

-  wa2Ap  +  Af/n  +  £  (ri+i-ri) [ r^ATi  +  ri+iATi+i]Ect/(l-2v)  (18) 

1=  1 

Now  we  can  form  a  system  of  2n+3  equations  for  solving  2n+3  unknowns,  (Aeg)^, 
(Aer)i,  at  i  -  1 ,2, . , .  ,n,n+l  and  A ez.  Equations  (16),  (17),  and  (18)  are 
taken  as  the  first  and  last  two  equations,  respectively,  and  the  other  2n 
equations  are  set  up  at  i  -  1,2,... ,n  using  equations  (13)  and  (14).  The 
final  system  is  an  unsymmetric  matrix  of  arrow  type  with  the  nonzero  terms 
appearing  in  the  last  row  and  column  and  others  clustered  about  the  main 
diagonal,  two  below  and  one  above. 

INCREMENTAL  LOADING  -  FIXED  VS.  SCALED 

When  the  total  applied  pressure  p  or  temperature  Tj.(i«l  to  n+1)  is  given, 
it  is  natural  to  divide  the  loading  path  into  m  equal  fixed  Increments  such  as 
Ap  -  (p-p*)/m,  ATi  -  (Ti~Ti*)/m.  These  fixed  increments  need  not  be  equal  for 
all  steps  and  any  sequence  of  m  increments  can  be  supplied  by  the  user.  A 
sequence  of  decreasing  load- increments  is  a  better  choice  than  that  of  equal 
Increments. 
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In  order  to  increase  the  efficiency  of  the  program,  an  adaptive  algorithm 
based  on  a  scaled  incremental-loading  approach®  has  been  implemented.  In 
each  step,  a  dummy  load-increment  such  as  Ap  is  applied  and  the  incremental 
results  Ao^  for  1  -  r,®,z  at  all  grids  are  determined.  For  all  grid  points  at 
which  o  ■  |  jojJ  |  <  oQ,  we  compute  the  scaler  “'s  by  the  formula 

o  .  i  {r  +  [T2  +  4| j Ao± { )2(a02-| | { \z)\l/2}/\ pOjJ |2  ,  (19) 

where 

r  *  ll°il|2  +  llAoil|2  -  ll°i  +  Aoi II 2  .  (20) 

and  ||°i||2,  ||Aoi||2*  li°i+  Aoil|2  are  computed  by 

IKH2  -  j  [(or-ae)2  +  (a0-a2)2  +  (az-°r>2J  •  (21) 

Let  X  be  the  minimum  of  the  a's.  Then  X  is  the  load-increment  factor  just 
sufficient  to  yield  one  additional  point.  A  sequence  of  X(J)  can  be 
determined  for  all  steps  j  -  i,2,...,m  and  the  updated  results  are 

p( j )  ■  p C j “1 )  +  X(J)Ap(J) 

oi(j)  -  o1(j-U  +  x(j>Ao1(J)  ,  etc.  (22) 

NUMERICAL  RESULTS  AND  DISCUSSIONS 

Consider  a  closed-end  tube  subjected  to  internal  pressure  p,  inner  and 
outer  temperature  T4  and  Tb»  The  numerical  results  were  based  on  the 
following  parameters:  b  ■  2",  a  ■  I",  n  ■  100,  v  ■  0.3,  E  ■  30xl06  psi,  do  * 
30xl03  psi,  «  -  0,  a  -  7.75xlO-6  in/in°F. 

According  to  Bland, ^  let  us  define 

6  -  EaN/t2(l-v)o0//3] 


’  me 4+ 


as  a  measure  of  the  effect  of  the  temperature  differences  in  the  stress 


system.  The  mean  temperature 


T\  ---  Trdr 


is  taken  as  zero  in  the  caicuiation  of  u/r.  Three  values  of  the  temperature 
stress  factor  8  are  chosen:  0  =  -1/2,  0,  and  1/2,  i.e.,  T\j-Ta  =  -72. 3°  F,  ' ) 3  F , 
and  72.3°F.  The  thermal  stresses  are  in  the  elastic  range.  In  the  presence 
of  these  temperature  gradients,  internal  pressure  p  is  applied  incrementally 
until  the  fully  plastic  state  is  reached.  The  internal  pressure  p  and  ins  1 
displacement  Ua  are  obtained  as  functions  of  elastic-plastic  interface  p  as 
shown  in  Figure  1.  The  effect  of  temperature  gradient  on  these  relations  is 
clearly  shown  in  the  figure.  In  order  to  compare  the  results  by  Bland,5  the 
state  of  stress  is  evaluated  at  p/a  =  1.73.  The  dimensionless  stresses 
(or  ,OQ,az)/aa  are  expressed  as  functions  of  r/a  as  shown  in  Figures  2  through 
4  for  three  cases  of  8  *  0,  -1/2,  1/2,  respectively. 

After  removing  the  temperature  gradients  and  internal  pressure,  the 
residual  stresses  were  obtained  for  all  three  cases.  The  states  of  residual 
stresses  for  the  first  two  cases  were  still  elastic  but  reverse  yielding 
occurred  when  unloading  the  last  case.  Five  scaled  incremental-loading  steps 
and  one  applied  incremental-loading  step  were  needed  to  unload  completely. 

The  residual  stresses  for  the  last  case  are  shown  in  Figure  5.  The  results  in 
Figures  1  through  5  are  roughly  the  same  as  those  of  Bland^  except  oz. 

This  difference  is  reasonable  because  Tresca's  yield  criterion  is  used  in 
reference  5. 
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As  a  laat  example  let  ue  consider  a  closed-end  tube  subjected  to  Inner 
temperature  Ta  only.  When  the  temperature  gradient  Is  of  sufficient 
magnitude,  yielding  will  first  expand  from  the  Inside.  When  Ta  Is  larger  than 
a  certain  limit  (238. 4*F),  plastic  zone  will  expand  from  both  the  Inside  and 
outside  surface  toward  the  Interior.  The  relation  between  the  Inside 
temperature  and  elastic-plastic  Interface  is  shown  in  Figure  6.  The  stresses 
in  a  closed-end  tube  subjected  to  temperature  gradient  (Ta  -  299° F,  Tb  •  0) 
are  shown  in  Figure  7.  The  dotb^I  lines  are  elastic-plastic  Interfaces. 
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gure  1.  Internal  pressure  and  displacement  as  functions  of  elastic-plastic  interface 
in  a  closed-end  tube,  (a)  Internal  pressure;  (b)  Internal  displacement. 
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Figure  3.  Straasaa  in  a  cloa*d-*nd  tub*  *ubj*et*d  to 
internal  pr***ur*  and  t«ap*ratur*  gradient 
(p/«  -  1.73;  0  •  -1/i). 
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Figure  4.  Stresses  la  s  closed-sad  tube  subjected  to 
laterael  pressure  end  teapersture  gradient 
(p/e  -  1.73}  6  -  1/2). 
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Figure  The  relation  between  internal  temperature 
and  elaatic-plaetic  interface,  •  0. 


Figure  7.  Stresses  in  •  closed-end  tube  subjected  to  teaperature  gradient. 
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